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Abstract 

Pions and protons production cross-sections are analyzed in proton-proton and 
proton-nucleus collisions at the RHIC energy at midrapidity. We employ the pQCD 
factorization scheme supplemented with the color-dipole formalism to investigate 
the Cronin effect. We calculate the broadening in the color-dipole approach for 
different centralities. Our main goal is to investigate, in a parameter-free man- 
ner within a unified framework, how much of the cronin effect for both pions and 
baryons comes from the transverse momentum broadening due to initial partons 
multi-scatterings. We conclude that final-state effects in pA collisions are impor- 
tant. Uncertainties in nuclear shadowing of various parton distributions and parton 
fragmentation functions are also discussed. 
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1 Introduction 



It is believed that p + A collisions provide a decisive testing ground to distin- 
guish between the initial- and final-state (plasma) effects in A + A collisions 
and can be used as a baseline for jet-quenching models. The Cronin effect 
[1] (which is generally associated with the ratio of p + A and p + p cross sec- 
tions, scaled by the number of collisions) in high-energy collisions has been 
the subject of renewed interest in recent years [2,3,4,5,6,7,8,9]. In order to pin 
down the role of parton energy loss effects in heavy ion collisions, a precise 
and firm understanding of the underlying dynamics of the Cronin and the 
shadowing effects in p + A collisions is indispensable. 
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There have been two very different approaches to explain the Cronin effect in 
p + A collisions: the initial-state effects [1,2] due to the broadening of the par- 
ton transverse momentum in the initial-state where the fragmentation of hard 
partons is assumed to occur outside the cold medium, and final-state effects 
[6] due to the recombination of soft and shower partons in the final-state. To 
understand the role of initial-state effects in the p + A reactions, one should 
also understand the broadening of transverse momentum of a projectile parton 
propagating and interacting with a nuclear medium. A first principle calcula- 
tion of the broadening of transverse momentum of partons is very complicated 
which involves apparently both soft and hard interactions. A common prac- 
tice involves fitting to a given experiment and then extrapolation to another 
reaction [3]. However, such an approach is less reliable since broadening is not 
universal and depends on kinematics and reaction. 



Here, we employ the pQCD factorization scheme supplemented with the color- 
dipole formalism to investigate the Cronin effect for identified hadrons. The 
broadening is calculated in a parameter-free way in the color-dipole approach. 
Our modest goal is to investigate the role of initial-state effects in the observed 
Cronin ratio for both protons and pions in p + A collisions and to learn how 
much of the effect comes from the broadening due to initial partons multi- 
scatterings. One of the uncertainties in such a calculation from the outset 
is due to our lack of knowledge about the baryon fragmentation functions. 
We investigate the implication of the uncertainties among various nuclear 
parton distributions and parton fragmentation functions and the role of parton 
primordial transverse momentum on the Cronin effect. 



One of the stunning experimental observations in both p + A and A + A 
collisions has been the much larger magnitude of the Cronin effect on baryons 
compared to mesons at the intermediate transverse momentum, the so-called 
baryon/meson anomaly in high energy nuclear collisions [10,11]. There have 
been different approaches to understand this phenomenon [6,12]. But a clear 
explanation is still lacking. Here, we compute the protons/pions ratio in p + p 
and p + A collisions. We also address if the underlying production mechanism 
of baryons and mesons are the same. 



This paper is organized as follows: In Sec 2, we discuss hadron production in 
p + p collisions. In Sec. 3, we introduce the main elements of our computation 
hadron spectra in p + A collisions. In particular, we discuss how to compute 
broadening in the color-dipole formalism. Sec. 4 summarizes our work. 
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2 Particle production in p + p collisions 



High-pT inclusive particle production in hadron scattering may be described 
by collinear pQCD factorization. However, in the derivation of the collinear 
factorization theorems, one applies to the hard scattering the approximation 
that the transverse momentum of the incoming parton can be neglected with 
respect to the transverse momentum generated in the scattering, and one also 
neglects the transverse momentum generated in the fragmentation. When the 
transverse momentum of the incoming partons is in the order of hard scale Q, 
the errors in the collinear approximations can be compensated by a correct 
treatment of higher order corrections to the hard scattering. However, for low 
transverse momentum compared with Q, in collinear factorization there is no 
precise compensation from higher order correction and one should not neglect 
the transverse momentum of the incoming partons. 

Moreover, it is well known that the curvature of the hadron spectrum can be 
corrected by considering intrinsic transverse momentum kx for the colliding 
partons. It was observed that significant parton intrinsic transverse momen- 
tum is needed for describing data of Drell-Yan dilepton production [13], direct 
photon production [14,15] and heavy quark production [16]. The origin of 
such a large primordial transverse momentum still needs a clear explanation. 
Nevertheless, some higher order pQCD corrections and soft gluon radiation 
corrections are effectively embodied in the primordial transverse momentum 
distribution. There are also some spin-dependent effects measured by experi- 
ments, such as single spin asymmetry, which cannot be explained by collinear 
factorization in all orders of QCD. The only way to produce these non-zero 
effects is with the presence of parton intrinsic transverse momentum [17] or 
by considering the higher-twist contributions. 

There exists many ways of incorporating the parton intrinsic transverse mo- 
mentum phenomenologically, and we choose for simplicity a kx smearing of 
the cross section to approximate this effect via unintegrated parton distribu- 
tions. This pQCD-improved parton model based on factorization should not 
be mixed with the collinear factorization. Notice that we do not aim at com- 
paring different types of factorizations to claim which one is better or not. 
Each of them can be applied to some specific processes and specific kinemat- 
ical region. Here, we prefer to use a kind of /c^-factorization which is flexible 
enough in the case of p + A collisions. 

By assuming the validation of the factorization in high-p^ particle production, 
and considering the intrinsic transverse momenta of the initial partons, the 
differential cross section in p + p collision can be written as [18] , 
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da pp ^ h+x r 
, , 2 = J2 I dx i dx j d 2 k iT d 2 k jT fi/ p (x u Q 2 )G p (k iT ,Q 2 ) 

y PT ijnl 

x f j/p ( Xj , Q 2 )g p (k ]T , Q 2 )K - nl) D ^ z ^\ (i) 

at Kz n 

where the K factor which is in general y/s and scale dependent accounts 
for the contribution of the NLO corrections (for the Cronin ratio which is the 
main subject of this paper, the K factor drops out and is not more important). 
fi/ P (xi, Q 2 ) is the parton distribution functions (PDF) of the colliding protons, 
which depend on the light-cone momentum fractions Xi and the hard scale Q. 
The function D h / n (z n ,Q 2 ) is the fragmentation function (FF) of parton n to 
the final hadron h with a momentum fraction z n . The differential cross section 
%(u — ► nl) of the hard process a{ki) + b(kj) — > c(k n ) + d(ki) can be calculated 
perturbatively(in the running coupling a s (Q 2 ) with scale parameter A equal 
to pion mass), in terms of the following Mandelstam variables 



k 2 k 2 

s = (h + kj) 2 = XiXjS + lT jT - 2k iT ■ k jT , (2) 

Ob iOb j S 

t = (ki - k n f = -l( XiV ^ PT e-y + -^prey - 2k tT ■ p T ), (3) 

Z n Xiy S 

1 k 2 
u = (kj - k n f = {xjyfe>re v + -^Pre^ - 2k jT ■ p T ), (4) 

Z n Xj y S 

and y^i is the center of mass energy of the incoming p + p system, y is the 
rapidity of the produced hadron defined as y — ^\n(E+p L )/(E—p L ). To avoid 
the divergence of the partonic cross-sections when one of the Mandelstam 
variables approaches to zero, we apply following replacement 

s^s + 2/i 2 , t^t — fi 2 , u-^u — fi 2 , (5) 

with fi = 0.8 GeV. This replacement will not change the partonic cross-sections 
at large s, t and u. Another advantage of this replacement is that the relation 
s + t + -u = 0is always satisfied. 

The function Q p (k a T) describes the distribution of intrinsic transverse mo- 
menta carried by partons. In phenomenological applications, a Gaussian form 
for Q p {k T ) is usually used: 

exp(-/4/(/4)) 
G ^ kT) = ^kj) ' (6) 

where (/c^) is the square of the 2-dimensional RMS width of the £>r distribution 
for one parton, which is related to the square of the 2-dimensional average of 
the absolute value of kr of one parton through (kj) = 4(/c T ) 2 /7r. For simplicity, 
we assume (k^) to be independent of Q 2 . 
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Fig. 1. Cross section of p + p -> vr° + X at RHIC energy <fs = 200 GeV for AKK08 
and Kretzer fragmentation functions. The error bars are the total uncertainties. 
Data are from the PHENIX collaboration [23]. 
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Fig. 2. Cross section of p + p — > tt + (tt^) + X and p + p — > + X at RHIC 
energy = 200 GeV. For theory curves we use the parameter set with AKK08 
fragmentation function, see the text. The error bars are the total statistical and 
systematic uncertainties. Data are from the STAR collaboration [10] . 
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Fig. 3. Ratio of (p + p)/(ir + + ir~) at RHIC energy y/s = 200 GeV for midrapidity 
with AKK08 FF. Data are for ratio of p/ir~ and p/ir + at midrapidity (| y |< 0.5). 
The error bars are the total statistical and systematic uncertainties. Data are taken 
from Ref. [10]. 



With the factorization form shown in Eq. (1), we calculate the cross section 
of p + p — > 7T° + X at RHIC for \/s = 200 GeV. For the parton distribution 
fi/ p (x,Q 2 ) used in our calculations we adopt the newest MRST PDF [19]. 
For the FF Dh/ n (z n ,Q 2 ), we first apply the parametrization by Kretzer [20]. 
However, Kretzer FFs do not give parton FFs to protons. Therefore, we also 
consider AKK08 parametrization [21] of FFs which give parton FFs to pions 
and (p+p). In both cases of calculations with Kretzer and AKK08 FFs, we 
will take (fcf.) = 2 GeV 2 and Q = Pt for the scale Q of the hard process. 
Notice also that different value of (k^) has been used in different approaches 
[9,22]. We take the if -factor K = 2.5, 1.5 for the cases of Kretzer and AKK08 
parametrization of FFs, respectively, which gives a good approximation of the 
higher order contribution in the pt region of interest. With the settings men- 
tioned above, we are able to reproduce the pions (tt°, n + + n~) productions 
data from the STAR and PHENIX collaborations [23,10] for p + p collisions. 
The discrepancy between theory and the data is within 40%, see Figs. (1,2). 
The setting with AKK08 FF also gives good predictions for protons produc- 
tions in p + p collisions in accordance with data from the STAR collaboration 
[10], see Fig. (2). 



In Fig. (3), we show experimental data from the STAR collaboration for ratio 
ofp/7r~ andp/7r + at midrapidity. The ratio of (p+p) / (tt + + tt~) at midrapidity 
obtained with AKK08 FF, seems to be consistent with the experimental data 
forp/7r~ and p/n + . 
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3 p + A collisions and Cronin effect 



Multiple interactions of projectile partons in the target may proceed coher- 
ently or incoherently. In the former case the multiple interaction amplitude is 
a convolution of single scattering amplitudes, and in the latter case one should 
convolute differential cross sections, rather than amplitudes. The underlying 
mechanisms of the multiple particle interactions and particle productions are 
controlled by the coherence length [24] , 

k ^ (7) 
m N p T 

where px is the transverse momentum of the fragmented hadron at midrapid- 
ity, and mjy is the nucleon mass. For pion production, the average momentum 
fraction (z) in the FFs is about 0.4 — 0.6 in the range of 2 < p T (GeV) < 8. For 
a coherence length which is shorter than the typical internucleon separation 
lc ^ Ra (where R A denotes the nuclear radius), the projectile interacts inco- 
herently. At the RHIC energy, y/s = 200 GeV, within intermediate pt we are 
almost in the transition regimes between the short- and long-coherence limit 
for more central collisions, and at higher p? we are in the short-coherence 
length (SCL) limit. 

Here, we resort to a simple scheme assuming that the pQCD factorization is 
valid at the RHIC energy at the midrapidity. We assume that the high-p^ 
hadrons are mainly originated from projectile's partons whose transverse mo- 
mentum is broadened by partons multi-scattering via gluons exchange. The 
broadening is computed in a parameter-free manner in the color dipole ap- 
proach. The single inclusive particle cross section in minimum-biased p(d) +A 
collisions can be written as 



j a p(d)A^h+X 



J2 J dxidxjd 2 k iT d 2 kj T T A (b) fi/ p {d){xi,Q 2 )Q P {k iT) b,Q 2 ) 



dyd 2 p T d 2 b ijnl 

x f j/A ( Xj , b, Q 2 )G p (k jT , Q 2 )K - nl) Dh/nM2 \ (8) 

at 7rz n 

where T A (b) is the nuclear thickness function normalized to / d 2 bTA{b) = A. 
We will use the Woods-Saxon nuclear profile for 7^(6). fj/A(xj, b, Q 2 ) denotes 
the parton distribution function (Npdf) in the target nuclei (with atomic num- 
ber A and charge number Z) which can be parametrized in a factorizable form 



\Z ( Z\ 

f j/A (xj, b, Q 2 ) = R j/A (x, b, Q 2 ) -f j/p (x, Q 2 ) + (l - — J f j/n (x, Q 2 



, (9) 



in terms of parton distribution in a nucleon fj/ n (x, Q ) and the nuclear mod- 
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ification function factor Rj/A{x,b,Q 2 ). We will show the results of our calcu- 
lation using three different Npdfs: EKS [25] and HKN [26] which are impact- 
parameter independent but the Q 2 -scale dependent via DGLAP evolution 
equation and HIJING-new [27] which is impact-parameter dependent. 

Initial/final state broadening of the projectile/ejectile partons is effectively 
taken into account via a modification of the primordial transverse momentum 
distribution, 

g p (k lT , b, q 2 ) = dat=q{q f^ qX \ x , b), (io) 

d 2 k iT 

where the transverse momentum distribution of partons after propagation 
through nuclear matter of thickness T^{b) is subjected to the broadening and 
computed in terms of the propagation of a qq color dipole through the target 
nucleus [28] 



da i=q (qA -> qX) 



d 2 p T " ' (2ny 

(11) 

where Qf n (fi,f2) is the density matrix which describes the impact parameter 
distribution of the quark in the incident hadron, 
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where (kj) denotes the mean value of the parton primordial transverse mo- 
mentum squared. For a projectile gluon (i = g) we replace the qq dipole with 
a gg one which can be written again in terms of qq dipole via a Casimir factor 
(j gg = \<Jqq- The appearance of the dipole cross section in Eq. (11), is the result 
of a product of the amplitude and the time-conjugated one, which describe the 
quarks with different impact parameters. Clearly, the object participating in 
the scattering is not a qq dipole but rather a single colored parton. To simplify 
the calculations we assumed that the initial and final partons are the same, 
so the total nuclear thickness T^ip) contributes to broadening. 

We take for the dipole cross-section in Eq. ( 11) for the case of a projectile 
quark (i = q), the popular saturation model of Golec-Biernat and Wusthoff 
[29]: a qq -(x,r) = a (1 — exp (— r 2 / B%)) where the parameters, fitted to DIS 
HERA data at small x, are given by a = 23.03 mb, R = 0.4fm x (x/x ) 0144 , 
where x = 3.04 x 10~ 4 . This parametrization gives a good description of DIS 
data at x < 0.01 [29]. 

The broadening in this scheme depends on the impact parameter b, the trans- 
verse momentum of incident parton k T , and also on the target Bjorken x. Note 
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Fig. 4. The mini-biased Cronin ratio R for ir° production at RHIC energy yfs = 200 
GeV. All curves are at midrapidity. In the upper panel, we show RdAu with AKK08 
FF for different nuclear PDF. In the lower panel, we show RdAu with HIJING-new 
nuclear PDF for different FFs. The normalization uncertainty on the p+p reference 
of 9.7% is not included here. The error bars are the total statistical and systematic 
uncertainties. Data are taken from Ref. [23]. 

that in Eq. (8), the summation over index i,j,n,l for quarks and gluons are 
different, not only because the corresponding pdfs and FFs are different, but 
also because the broadening of the projectile quarks and gluons are different. 



The i^-factor in Eq. (8) and the parton primordial transverse momentum 
squared (k^) in Eqs. (8,12) are taken the same value fixed in p + p collisions. 
Therefore, all the phenomenological parameters in the master equation (8) are 
already fixed in reactions different from p(d) + A collisions. 

The nuclear modification factor R p a defined as ratio of p + A to p + p cross 
sections normalized to the average number of binary nucleon collisions, 



da pA^h+X 

p _ dyd?p T / 1C >\ 

V* binary/ d yd 2 p T 



where (Nonary) is the average number of geometrical binary collisions which 
is calculated according to the Glauber model [30] for different centrality. 

In Fig. (4) (upper panel), we show mini-biased RdAu for n° production at the 
RHIC energy for deuteron-Gold collisions with AKK08 FF for various Npdfs. 
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Fig. 5. The Cronin ratio R for ir° production at RHIC energy y/s = 200 GeV at dif- 
ferent centralities. All curves are at midrapidity with AKK08 FF and HIJING-new 
nuclear PDF. Data are taken from Ref. [23]. 



It is obvious that the uncertainty among various Npdfs at small px in the 
shadowing region leads to rather sizable different Cronin ratios. The position 
of the RdAu peak in pt does not seem to vary that much among various Npdfs. 

In Fig. (4) (lower panel), we show RdAu for n° production with HIJING-new 
Npdf for different FFs. As one may expect both considered fragmentation 
functions (Kretzer and AKK08 FFs) give very similar results for RdAu, since 
we assumed no medium modification for FFs. Notice also that the energy loss 
is less important at the RHIC energy for d + Au collisions, although it might 
be important at lower energies at SPS [5,31]. 

One should note that the nuclear interactions that lead to shadowing are also 
the source of parton momentum broadening. We incorporate the shadowing 
effect and the other nuclear medium modification of the partons through the 
Npdfs Eq. (9). In principle, at high-energy the qq dipole cross-section used for 
obtaining the broadening is also subjected to the multi-Pomeron fusion effect 
in the presence of nuclear medium. However, we have already incorporated 
the shadowing effect in the Npdf and including the shadowing effect into the 
dipole cross-section may lead to double counting. One should also note that the 
onset of gluon shadowing and its magnitude are still debatable, see Ref. [32] 
and references therein. 



In Fig. (5), we show the Cronin ratio at different centralities at midrapidity. 
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Fig. 6. The mini-biased Cronin ratio R for charged pions 7r + + 7r~ and baryons p + p 
productions at RHIC energy y/s = 200 GeV at midrapidity with AKK08 FF for 
various Npdfs. The experimental data are for | y |< 0.5. The shaded bands around 
1 (of the order of 17%) corresponds to the error due to uncertainties in estimating 
the number of binary collisions in minimum bias. Data are taken from Ref. [10]. 
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Fig. 7. The mini-biased ratio of {p + jj)/(vt + + ir~) for d + Au collisions at RHIC 
energy = 200 GeV for midrapidity with AKK08 FF and various Npdfs. Data 
are for ratio of p/ir~ and p/ir + at midrapidity (| y |< 0.5) for minimum bias d + Au 
collisions. The error bars are the total statistical and systematic uncertainties. Data 
are taken from Ref. [10]. 

We used for all curves AKK08 FF and HIJING-new Npdf. It is obvious that 
generally the centrality dependence is rather weak in accordance with the 
PHENIX data [23]. 
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Fig. 8. The mini-biased Cronin ratio for pions and protons at midrapidity for the 
RHIC energy with AKK08 FF and HIJING-new Npdf. 



In Fig. (6), we show the mini-biased RdAu at midrapidity for charged pions 
7r + + 7r~ and baryons p + p at midrapidity with AKK08 FF for various Npdfs. 
It is seen from Figs. (4,5,6) that our results obtained with HIJING-new Npdf 
agree with the Cronin data for both n° and charged pions ix + + 7r~. However, 
the same setting does not seem to be in a good agreement with the Cronin 
data for baryons in p + A collisions although it gives a good description of 
baryons cross-section in p + p collisions. Taking the experimental data at face 
value, the discrepancy between the theory and experimental data seems to 
persist up to px = 6.5 GeV. A number of studies have recently found that 
recombination is more important than fragmentation at small and moderate 
Pt at midrapidity in heavy-ion collisions at the RHIC energy ^Js = 200 GeV 
[6,33]. Therefore, a precise experimental measurement of the Cronin ratio for 
baryons at higher px > 7 GeV is essential in order to reveal the underlying 
mechanism of hadron production in the cold nuclear matter. This deviation is 
also seen from Fig. (7) where similar to Fig. (3), we plot (p + p)/(7r + + 7r~) 
ratio in d + Au collisions at midrapidity for various Npdfs. From Fig. (7), it 
is obvious that baryons/pions ratio in p + A collisions is not sensitive to the 
shadowing effects in contrast to R p a- To conclude, these discrepancies indicate 
that not entire Cronin effect in p + A collisions comes from the broadening 
due to initial partons multi-scatterings, and final-state effects which are not 
included here, are also important. 

In Fig. (8), we compare the mini-biased Cronin ratio for pions and protons at 
midrapidity for the RHIC energy with AKK08 FF and HIJING-new Npdf. It 
is obvious that the Cronin ratios for baryons and mesons coming from initial- 
state effects are not significantly different. 
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4 Summary and final remarks 



In this paper, we studied pions and baryons productions in p + p and d + 
Au collisions at the RHIC energy. Having a successful description of pions 
and protons productions in p + p collisions, we showed that the observed 
Cronin ratio for pions in d+Au collisions can be fairly described by transverse 
momentum broadening due to initial partons multi-scatterings. But the same 
mechanism seems to underestimate the observed Cronin ratio for protons in 
d+Au collisions. We stress that all phenomenological parameters in our model 
are fixed in reactions different from d + Au collisions. This discrepancy might 
indicate that initial state-effects (assuming that the fragmentation of hard 
partons occurs outside of medium) might not be totally accountable for the 
observed Cronin effect for baryons in p + A collisions and the separation of 
partons into two non-interacting components, soft and hard, might be an 
oversimplification. It is also possible that the baryons productions mechanism 
in cold nuclear medium to be different from pions productions. 

A similar approach to this paper, was taken by Kopeliovich et al. in Ref. [5] 
to calculate the mini-biased Cronin ratio for tt° . In this paper, we extended 
their study by calculating the Cronin effect for the charged pions and protons 
at different centralities. One of the major difference between our study with 
Ref. [5] is the way that the shadowing is included. As we already mentioned at 
moderate pt at the energy of RHIC, we are in the transition regime between 
the short- and long- coherence limits. Here, we assumed that the pQCD fac- 
torization is still valid in the transition region, and we included the shadowing 
effect in the conventional way via the nuclear parton distribution. We also 
studied the effect of various available nuclear parton distributions. However, 
in Ref. [5] no shadowing was included in Npdf. While the Cronin ratio was 
obtained by a linear interpolation between two Cronin ratios obtained in two 
different schemes of the short- and long-coherence limits. We also tried their 
prescription in order to investigate the pions and protons productions in d + A 
collisions. We found similar results to those obtained in our approach which 
leads us again to the same conclusion that baryons and pions productions 
mechanism in high-energy p + A collisions may be different and final-state 
effects are important [6,33]. Nevertheless, given rather large experimental un- 
certainties further studies are needed in order to make a final verdict. 

Notice that our prescription is not reliable at forward rapidities. The rea- 
son is due to the importance of the valence quarks contribution which are 
not incorporated in the color dipole picture. It might be tempting to assume 
that forward rapidity region should be valid domain of our scheme since the 
Bjorken X2 of target is small, and the color-dipole approach which is based on 
Pomeron-exchanges should be at work. However, one should note that at this 
region x± — > 1 and as a consequence the energy conservation put a restricted 
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constraint on the particle productions and therefore valence quarks become 
very important [8]. 
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